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Number counts of massive high-redshift clusters provide a window to study primordial non- 
Gaussianity. The current quality of data, however, forces the statistical analysis to probe a region 
of parameter space - the extreme tail of the mass function - which is neither accessible in any of 
the currently available theoretical prescriptions for calculating the mass function, nor calibrated in 
Ai'-body simulations. In this work we present a new analytical prescription for calculating a "re- 
summed" non-Gaussian halo mass function, which is constructed to remain stable in the extreme 
tail. We show that the prescription works well in the parameter regime that has been currently 
explored in simulations. We then use Fisher matrix techniques to compare our prescription with 
an extrapolated fit to A''-body simulations, which has recently been used to obtain constraints from 
data collected by the South Pole Telecope. We show that for the current data, both prescriptions 
would lead to statistically consistent constraints. As the data improve, however, there is a possibil- 
ity of introducing a statistically significant bias in the constraints due to the choice of prescription, 
especially if non-Gaussianity is scale dependent and becomes relatively large on cluster scales. It 
would then be necessary to test the accuracy of the prescriptions in A-body simulations that can 
probe clusters with high masses and redshifts in the presence of large non-Gaussianity. 



I. INTRODUCTION 

One of the key aims in cosmology is to determine hovif much non-Gaussianity there was in the primordial curvature 
fluctuations. Canonical single field inflation with a smooth potential predicts that primordial non-Gaussianity (NG) 
should be extremely small . Detecting an appreciable amount of primordial NG could therefore rule out these 
models. Measurements of the cosmic microwave background (CMB) [4] and the scale-dependent halo bias [5] provide 
good constraints on primordial NG on scales of about k ~ 0.05/i~^Mpc. Although the Planck experiment will improve 
the range of scales constrained the primary CMB will alway be limited by foregrounds to < 0.2/i~^Mpc. Since 

it is possible for primordial NG to be quite strongly scale dependent [1, , it is important to constrain it on as wide 
a range of scales as possible, with complementar y p robes when possible. Number counts of galaxy clusters provide 
such a probe on smaller scales (see also Refs. Il0l-[T3). There have been several recent attempts at using massive 
high-redshift clusters to constrain primordial NG llill^. Wilhamson et al (Wll) [13 analysed a subset of the 
clusters detected by the South Pole Telescope (SPT), using a likelihood analysis which carefully accounts for issues 
such as the survey selection function. They find a posterior probability distribution for the standard non-Gaussianity 
parameter /nl, and their quoted result is /nl = 20 ± 450 at 68% confidence. 

A key ingredient in any such analysis is the chosen prescription for the non-Gaussian halo mass function. Given a 
mass function, one can construct a likelihood for the data by computing the expected number of clusters in a given 
mass and redshift range, as an integral of the mass function. The Wll analysis uses the prescription of Dalai et al. 
(DOS) [m, which is essentially a fit to A^-body simulations. As we will see however, the Wll likelihood function 
probes a region of parameter space (the extreme tail of the mass function) which has not been calibrated by DOS (nor 
indeed, by anyone else). Additionally, all the other currently available prescriptions [tI. [l9l - l22| formally break down 
in this extreme region. 
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It is therefore interesting to ask whether one can obtain any analytical understanding of the mass function in 
its extreme tail, and whether the DOS prescription might be introducing a bias in the analysis. The main result 
of this work is a new prescription for the non-Gaussian mass function, which involves a resummation of an infinite 
perturbative series. The resulting expression compares well with the results of iV-body simulations. More importantly, 
it does not formally break down in the region of parameter space that needs to be probed by the Wll data analysis, 
and can therefore be used to compare against the DOS prescription. Using a Fisher analysis, we find that for the 
current quality of data the two prescriptions give statistically identical results. With better quality data however, we 
show that there can be a statistically significant bias between the two methods, especially if /nl on small scales is 
large. It would then be necessary to accurately calibrate the tail of the non-Gaussian mass function in simulations. 

The paper is organised as follows : Section |lT] motivates and introduces our new "resummed" prescription for the 
mass function, comparing it with other prescriptions, both in regimes which have been tested by simulations and 
in the extreme tail which has not. In section IIIII we use a Fisher analysis to compare the error bars on /nl from 
the DOS and resummed prescriptions, and analyse the level of bias between the two methods. We end with a brief 
discussion in section |lVl Technical details of calculations have been relegated to the Appendix. Unless otherwise 
specified, we assume a fiat ACDM cosmology with parameter values compatible with WMAP7 [3] : = 0.703 
with Hq = 100ftkms~^Mpc^^ the Hubble constant, total matter density flmh^ = 0.134, baryonic matter density 
ilbh^ = 0.0227, scalar spectral index Ug — 0.966, and erg — 0.S09. Throughout we use the transfer function of Bardeen 
et al. ^2^ , with a baryonic correction as prescribed in Ref. [24j . 



II. NUMBER COUNTS IN THE EXTREME TAIL 

Within the paradigm of hierarchical structure formation, any physically acceptable mass functior0 must be compatible 
with the fact that massive objects tend to form late in the cosmic evolution. Mass functions derived assuming Gaussian 
initial conditions 25-27]easily fit the bill by virtue of being monotonically decreasing at large masses and redshifts. 
The situation with non-Gaussian initial conditions is not as clean however. There are several prescriptions in the 
literature for theoretically calculating the non-Gaussian mass function 

0, El-il- All of these rely on a perturbative 
treatment of non-Gaussianity and assume a condition which is at least as strong as \ev\ < 1 (see D'Amico et al. (Dll) 
[2H for a detailed comparison), where e ^ /nlV^ is a small parameter controlling the non-Gaussianity {A ^ 10~^ 
being the power spectrum normalisation) and 

iy{M,z) = — — , (1) 

with the threshold Sc — 1.686 (for spherical collapse, but see below), the linear growth rate D{z), and the variance 
aj^ = (Sj^) oi the linearly extrapolated density field smoothed on scale R = (3M/47rp)^/'^. While several of these mass 
functions have been tested in A'^-body simulations and do reasonably well in their regime of validity, when extrapolated 
to \eiy\ > 1 they lead to unphysical results (e.g. the Matarrese et al. (MVJ) [l^ result becomes imaginary for /nl > 
when ev ~ 3). This would be an academic issue, were it not for the fact that the current quality of data requires 
predictions of number counts at combinations of {f-f^i^, M, z) values which lie squarely in the |e:/| > 1 corner [16]. To 
see this, note that e.g. the Wll analysis proceeds by calculating the joint probability, for a given /nl, of observing 
the ensemble of clusters in their sample, and of observing nothing else in the parameter range they explore. The large 
errors on /nl (~ 450) indicate that their likelihood is non- negligible at say /nl ^ 750. Fig. [T] shows eiz/ as a function 
of mass at redshift z = 1.5 for some representative /nl values, where 

ei^{Sl,)/al,, (2) 

with £1 ~ 3 X 10^'^/nl for a WMAP7 cosmology, and we identify ei with the small parameter e mentioned above. 
(Throughout this paper we will use the local model of NG to compute ei, see e.g. Dll for details.) We see that for 
/nl — 750 at this redshift, ev ^ 1 at M ^ 5 x lO^'^ft^^M©, which is well inside the region explored by Wll. 

In contrast with the theoretical approaches, the DOS prescription constructs a non-Gaussian mass function by 
convolving a Gaussian mass function with a Gaussian probability distribution whose mean and variance depend on 
/nl- As it turns out, this prescription results in the only viable mass function currently on the market which remains 
stable at large {fyh, M, z). Since the Wll analysis uses this prescriptiorlfl, it at least does not suffer from the 



^ By mass function we mean f^j^y (dV/ dz)(dn / dM); where dn is the comoving number density of halos with masses in (M, M + dM), f^^y 

is the fraction of sky observed and dV is the volume element with dV/dz = 4:TrH(z)~^ [J^ dz' H{z')~^'\^ . 
^ Presumably Wll used the exact numbers which DOS quote for the parameters of their fit; this detail is not mentioned in Wll. 
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FIG. 1: The quantity eiv, where ei and v were defined in Eqns. ((2| and IH) respectively. The curves are for redshift z 
and three difi'erent values of /nl, all probed by the Wll likelihood analysis. 



1.5 



breakdown which the MVJ, Lo Verde et al. (LMSV) 't'I or LoVerde & Smith (LS) [131 mass functions would encounter. 
There is however some cause for concern. Firstly, the actual parameter values DOS quote are a numerical fit to 
simulations which had considerably large Poisson errors and scatter at large M and z (see Fig. 6 of DOS). Subsequent 
simulations by Pillepich et al. (PPH) [2^ indicate that this specific fit does not work very well at the low mass end 
either. The main problem however is that this fit applies in a mass and redshift regime where |ez/| is at most ~ 0.5, 
and there is no reason to expect it to work accurately when ev > \. One would therefore like to have some analytical 
understanding of this extreme tail of the mass function. 



A. A New Stable Non-Gaussian Mass Function 



In this work we present a new analytical non-Gaussian mass function which is stable in the \ev\ > 1 regime. As we 
show in the Appendix, under some technical assumptions regarding the form of primordial non-Gaussianity, which 
are inspired by a known perturbative hierarchy in the local model, the perturbative series appearing in the analysis 
can be resummed. The excursion set formalism then gives a mass function which is valid at arbitrarily large v for 
a given positive /nl (which can also be large, see below). The formal excursion set result needs to be modified in 
order to predict the correct Gaussian mass function when /nl 0. Assuming that the Gaussian mass function is 
well described by the Tinker et al. form |2^, the "resummed" mass function we prescribe is 



dn 

dM 



Rcsum 



dn 
dM 



Tinker 



where the ratio 7?.Rosum is given by 



(3) 



Rcsum 



(M,z,/nl) = (l + eiu) 



-1/2 



exp 



-I- [eiv - (1 + eiv) ln(l -I- eiv)) 



(4) 



The expression is formally valid for {\-\-£iv) > 0, and receives a relative correction of order 0{e\{l +eiv)~^) due to a 
saddle point approximation (see the Appendix for details). This shows that although /nl cannot be made arbitrarily 
large for a given values of the order /nl ^ 1000 can be easily accomodated over the full (M, z) range, and the saddle 
point approximation becomes increasingly accurate for large (M, z). Finally, since the Tinker et al. mass function 
falls off like ~ g-o.837i5^/2CT ^ need to redefine the v that appears in 7?.Resum by replacing 



5c — !• V ^Tinker 



(5) 



in Eqn. ([T]) where ^Tinker = 0.S37. This ensures that the resulting mass function is still stable for arbitrarily large 
values of v. This correction is similar to the one first introduced by Grossi et al. [30|] and can be motivated by appealing 
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FIG. 2: Non-Gaussian ratios (dn/ciM)|NG/(rfn/dM)| Gauss as per different prescriptions, for values of /nl tested in the PPH 
A'^-body simulations These plots use the same cosmology as used by PPH in their simulations (the WMAP5 values in 
Table 2 of PPH). The curves correspond to resummed with q = gxinkcr ~ 0.837 (solid), LMSV [3| with q — 0.79 (dot-dashed), 
MVJ [T^ with q = 0.79 (dotted), and DOS ^ISl] (dashed) calculated using the Tinker et al. [2^ Gaussian mass function. The 
different sets of curves in each panel correspond to different redshifts. The lower panels show, from left to right, redshifts 
z — 1, 0.5 and 0. The upper panels show z = 1 and 0, with z — 1 corresponding to the lower set of curves in the left panel 
(/nl = —80) and vice- versa in the right panel (/nl = -1-80). We see that in this regime of parameter space, the resummed 
prescription (which has no free parameters once the Gaussian mass function is chosen) lies very close to LMSV with q = 0.79. 
These plots can be directly compared with Fig. 5 of PPH, who showed that LMSV with q = 0.79 fares very well compared to 
simulations. Consequently, so does the resummed prescription. 



to a stochasticity in the spherical collapse barrier [31| . We emphasize that this modification, and the specific value 
^Tinker, IS forced OH US by our choice of the Gaussian mass function - we have no free parameters to play with in our 
derivation. 

Fig. [5] shows the non-Gaussian ratios according to various prescriptions, for values of /nl tested in the PPH 
simulations. For better comparison with PPH, these plots assume the same cosmology that PPH used in their 
simulations (the WMAP5 values in their Table 2). We see that in the range plotted, the ratio 72.Rosum is very close to 
the LMSV ratio TZlmsv (which is plotted using a v corrected with q = O.TOj^. This can be understood analytically 
since in the regime eiy^ < 1, 7?.Rosum and 7?.lmsv become formally identical when truncated at order 0{ei^). The 
residual difference is due partly to higher order terms and partly to the different values of the g-correction. Fig. [2] 
can be directly compared with Fig. 5 of PPH, who showed that 7?,lmsv with q — 0.79 performs very well compared 
to simulations in this range. Consequently, so does our resummed ratio. 

For large negative /nl (or large v) such that ev < —1, our prescription breaks down. Nevertheless, there is a very 
simple way of extending the mass function to this regime. We notice that for |e^| <Si 1, it is approximately true that 
''^Rcsum(Af, z, — /nl) — "T^Rcsuml-W, z, /nl)~^- We therefore define the resummed ratio for negative /nl and arbitrary 
1/ by 

7?.Rcsum(Af, Z, -/nl) = CR-RcsumiM, Z, /nl))"^ ; /nL > . (6) 

Of course this is completely ad hoc, but we see from Fig. [2] that the prescription compares extremely well with 7?.lmsv 
at values of /nl < that were tested in the PPH simulations. It will be very interesting to see how the resummed 
ratio compares with simulations at high masses and redshifts, for both positive and negative /nl values. 



The MVJ and LMSV ratios in our language are respectively given by T^mvj = (l — i^^^) ^''^ e-i" (l — ^^iJ^ (l — §^i)) ^'^'1 
■^^LMSV = 1 + i^i'^'^ ( 1 V (3 — 2ei) ^ei j, where si = dlnsi/iilno-^. 
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FIG. 3: Non-Gaussian ratios (dn/ciM)|NG/(dn/dM)|Gauss as per difTerent prescriptions, for /nl ~ 500 and z = 1.5. 
Comparing with Fig. [T]we see that the curves now enter the regime ev > 1. The curves correspond to resummed with 
q = gxinkor = 0.837 (sohd), LMSV (dot-dashed) and MVJ (dotted) with q = 0.79, DOS (dashed) with the Tinker et al. 
Gaussian, and the "log-Edgeworth" ratio prescribed by LS [S^] (cross-dashed). The DOS prescription continues to 
systematicaUy predict a larger number of massive halos than the resummed. 



III. CONSTRAINTS ON /nl 

Fig. [5] also shows that for a given positive /nl, the DOS ratio 7?.Daiai is consistently higher than 7?.Rosum and therefore 
predicts more high mass halos. This trend remains true even in the regime ev > 1, as Fig. [3] shows. It is then worth 
asking how sensitive the analysis of clusters a la Wll is to this difference^. 

We use a Fisher analysis to make a "forecast" for the kind of data that SPT has already observed (specifically 
the subset used by Wll). We set our fiducial cosmology to the WMAP7 values given earlier, with /nl = 0. We 
marginalise over erg with a WMAP7 prior, but assume perfect knowledge of all other cosmological parameters. This is 
not expected to significantly affect our conclusions. We bin in redshift between 0.3 < z < 2, with a spacing Az — 0.05. 
The lower limit is the same used by Wll, and we have checked that the analysis is insensitive to increasing the upper 
limit. We assume that the survey sees all clusters in any given redshift bin above a certain constant mass threshold 
-A^iim, which is approximately true for a Sunyaev-Zeldovich survey such as SPT [33l |. We marginalise over A/iim with 
a lognormal prior of 30% in Miim, in keeping with typical mass uncertainties quoted by Wll. The fiducial value for 
Miim is chosen such that the total expected number of clusters with masses above and at redshifts z > 0.3 in 

the fiducial cosmology, is roughly the same as in the Wll analysis, which is 26. For /nl = 0, (i.e. - using the Tinker 
et al. mass function) this gives us a fiducial value — 8.5 x lQ^'^h~^ Mq for /sky = 0.06. 

We then assume a joint likelihood given by a product of independent Poisson probabilities for each redshift birlf], 
and construct the Fisher matrix for parameters Qa — (Afiim, cts, /nl) [35| 

where /ii is the expected number of clusters in the j*^ redshift bin, evaluated as the integral of {dn/ dM){dV / dz) over 
this bin, for masses above Mum. The results are in Table IH which shows the marginalised errors on /nl- We see 
that the DOS mass function results in a marginal error Of^-^ ~ 300, which can be compared with the error quoted 
by Wll which is 450. This shows that our analysis works; the fact that our marginal error for DOS is smaller, is 



^ We note that in contrast to Wll, earlier analyses Il5l |16(| showed a significant tension with standard ACDM. It would then seem to 
be more interesting to perform our analysis on the set of clusters used in these analyses rather than Wll. In fact, these analyses can 
be shown to have used a biased statistic, and removing this bias makes them consistent with the analysis of Wll |32| . We ultimately 
choose to study Wll because we find it easier to approximate the selection function for the Sunyaev-Zeldovich clusters measured by 
SPT. 

® This is similar to what is done by Wll, who use a product of independent Poisson probabilities over a number of bins in the space of 
redshift and detection significance which is roughly their mass proxy. As can be seen from the methods of Ref. [3^ , sample covariance 
should be negligible for this survey. 
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Marginal cr/j^^^, Wll equivalent 


Prescription 


/nl = 


/nl = 500 


DOS 


298 


463 


Resummed 


457 


599 


MVJ 


487 


288 


LMSV 


498 


1366 



TABLE I: The marginal Fisher error (Jf^j^ as per different prescriptions, for two different choices of fiducial /nl, for a survey 
that is approximately equivalent to the subset of the SPT clusters analysed in Wll (see text for details). The marginalisation 
is over ag with a WMAP7 prior and the threshold mass Mn^ with a 30% lognormal prior (see text). We see that the 
resummed and D08 prescriptions predict comparable errors in all cases, while the MVJ and LMSV predict significantly 
different values. See text for a discussion. 



(jjj^LJ Planck-like 




Conditional 


Marginal 


Prescription 


/nl =0 


/nl = 500 


/nl =0 


/nl = 500 


D08 


9.1 


9.2 


45.7 


77.9 


Resummed 


11.0 


11.0 


64.0 


80.5 


MVJ 


10.8 


6.6 


72.4 


41.7 


LMSV 


12.0 


16.1 


69.6 


310 



TABLE IL Fisher errors on /nl for a final Planck-like survey, with a 2% prior on and a fiducial threshold 
Afiim = 5 X 10^^h~^ Mq with a 10% lognormal prior and /sky = 0.8. The values under "conditional" assume perfect knowledge 
of all parameters except /nl, and represent what is possible if parameter degeneracies are broken using measurements of say, 
the clustering of clusters p^ffi^ . 



consistent with our simplifying assumptions regarding the cosmological parameters and selection function. The errors 
predicted by the resummed mass function are comparable to that predicted by the DOS prescription, indicating that 
using either prescription would lead to similar results. Since the errors are large however, one might wonder if the 
situation would change if the fiducial /nl were some large value consistent with the error. Table U shows that the 
predicted errors from the DOS and resummed mass functions are still comparable at a fiducial /nl — 500 (in this case 
we used a fiducial i\/iim = 1-05 x 10^^ Mq to get a total expected number close to 26). 

Table HIl shows the results of a Fisher forecast for a final Planck-like SZ survey [s^, assuming a 2% prior on erg, 
a threshold = 5 x 10^^h~^MQ with a 10% prior, and binning in redshift between 0.1 < z < 2 with a spacing 

Az — 0.05 and /sky — O.S. In this case we also display the conditional error (i.e. - assuming perfect knowledge of 
other parameters), as being representative of what can be achieved when parameter degeneracies are broken using 
measurements of say, the clustering of clusters |1( ML3| . We see that the predicted conditional and marginal errors for 
/nl are again comparable for the two mass function^. 

It appears therefore, that using the DOS fit at ei^ > 1, as opposed to a theoretically motivated mass function such 
as the resummed, does not introduce a significant effect in the error predicted for /nl from number counts of clusters. 
However, since the actual number of halos predicted at any given /nl is different for these two mass functions, we 
must also worry about a possible bias in the central value of /nl- Estimating this properly would require a full- hedged 
Monte Carlo analysis that accurately accounts for survey selection functions, which is work in progress. For now, 
we perform a cruder analysis. For parameters corresponding to the final Planck-like survey mentioned above, we 
assume a fiducial cosmology of the WMAP7 ACDM, but this time with a non-zero value of /nl = /nl* using the DOS 
prescription. We then "analyse data" using the resummed prescription, by constructing an approximate likelihood 



''' This might seem surprising given the difference between the magnitudes of these mass functions. In fact, it might appear from Fig. 
[3] that the MVJ prescription should be closer to DOS. Note however, that the Fisher matrix (0 depends not only on the integrated 
mass function, but also on its derivative with respect to /nl, leading to a complicated interplay which is difficult to predict simply by 
examining plots of the mass function or non-Gaussian ratio. 



7 




0( 



/nl 



FIG. 4: Testing for bias between DOS and resummed for a final Planck like survey : Posterior probability distribution 
P(/NL|/NLt) computed by marginalising the approximate likelihood ^ for the DOS "data" given the resummed "model", over 
Miim with a 10% lognormal prior. The curves are for /nl, ~ 100 (solid, black) and /nl, ~ 500 (dashed, red). The 
corresponding vertical lines mark the value of the respective /nl,- We see a statistically significant relative bias between the 
resummed and DOS prescriptions for large /nl,, while at smaller /nl, the prescriptions are statistically consistent with each 
other. 



jiven by 



'C(/NL,Mim) oc exp 



- E 



.Rcsum 



Dalai V 



.Rcsum 



In/i: 



Rcsum 



(8) 



where /i^*^'*^' — M?'^''^'(/nl,, Mim*) is the fiducial expectation value in the i^^ bin, which we treat as "data", and 
^Rcsum _ fif°^^'^{fNL, Mii^) IS the model expectation value computed using the resummed prescription at the 
(/nLj -^^lim) values being analysed. (We keep erg fixed in this exercise.) We marginalise this likelihood over Miim 
with a 10% lognormal prior (assuming a flat prior on /nl), to obtain a posterior probability distribution for /nl- Fig- 
m shows the results for /nl* = 100 and 500. We see that for /nl* = 100, the resummed posterior p(/nl|/nl* = 100) 
peaks close to /nl, and there is no statistically significant bias; /nl* lies less than one standard deviation away 
from the peak value. For /nl* = 500 on the other hand, the value /nl, lies far in the tail (more than 4 standard 
deviations away from the peak) of the distribution p(/nl|/nl* = 500), indicating a statistically significant bias. A 
similar calculation for the SPT equivalent parameter values used in Table [J shows that with the current data there 
would be no detectable difference between analyses based on the resummed or DOS prescriptions. We can conclude 
that systematic effects from the choice of non-Gaussian mass function could become important as the data improves, 
if non-Gaussianity is scale dependent and results in a large /nl on cluster scales. In this case it would be necessary 
to properly calibrate the non-Gaussian mass function in its extreme tail. 



IV. SUMMARY & CONCLUSIONS 



Number counts of clusters offer a very interesting probe of the extreme tail of the mass function of collapsed objects, 
and thereby the statistics of the primordial curvature fluctuations. In this work we developed a new prescription for 
calculating the modification of the mass function in the presence of primordial non-Gaussianity (Eqns. ^ and ([3])). 
The key features of this "resummed" prescription are that it is theoretically well motivated (as opposed to a fit to 
simulations like DOS ^ISJ) as well as stable in the extreme region |ez/| > 1 (unlike all the other currently available 
theoretical prescriptions 0, [13413 )• We also showed that this resummed prescription compares very well against the 
results of A^-body simulations in the literature (see Fig. [5]), despite having no free parameters left once its Gaussian 
limit is fixed as the Tinker et al. [29] mass function. 

While this might appear to be an academic issue, we showed that the current quality of data forces any likelihood 
analysis to enter the regime |ei^| > 1. It is then important to analyse how far one can trust the chosen prescription 
for the non-Gaussian mass function. The specific example we dealt with was the Wll [l3| analysis of a subset of the 
current SPT clusters, which used the DOS prescription to compute the likelihood. We used Fisher matrix techniques 
to analyse the consequences of changing the prescription from DOB to resummed, and found that for the current 
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data, both prescriptions are expected to give approximately identical results. As the data improve however, there is a 
possibility of introducing a statistically significant bias in the analysis, due to the choice of prescription. It will then 
become important to evaluate which is the better method, which will require improved A^-body simulations. This 
would be especially true if non-Gaussianity is scale dependent and results in a large value of /nl on cluster scales 
1^] . One would then need to calibrate the mass function using A^-body simulations which may have to be specially 
tailored (with larger particle size, say) to probe clusters with high masses, at large redshifts and in the presence of 
large /nl. 
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Appendix A: The Resummed Mass Function 



In this Appendix we sketch the derivation of the resummed mass function presented in the text, describing all 
the approximations that enter. The non-Gaussian halo mass function derived from the excursion set formalism for 
spherical collapse is 



dn 
dM 



— J— f 



diner 



dlnM 



NG.cxc 

where the multiplicity /ng.oxc for large masses is given by (see Dll [21]) 



(Al) 



/NG,oxc(i^,ei,e2, ■■■) ^ \l exp 



(-1)" 



.r!=3 



(A2) 



where z/(M, z) was defined in Eqn. ([T]) and £„_2 = {^m)/'^m normalised connected moments of the linearly 

extrapolated, smoothed density field. This expression assumes that the functions ej are all constant with scale in 
the regime of interest, which is reasonable at least in the local model (see e.g. Fig.l of Dll). If /nl is the only 
NG parameter allowed, then these are perturbatively ordered : ej ~ (/nl'v/^)"' , with A the normalisation of the 
power spectrum. The ellipsis in Eqn. (jA2p indicates terms arising from "unequal-time" correlations as discussed in 

Refs. [13, [2l|. In the large mass (or large v) limit, these will be of the form ~ e~i'^^ 0{ei^). Experience with the 
calculation of Dll indicates that the action of the exponential derivative on all these terms will result in a single 
common exponential prefactor multiplying a series expansion in ev. When \ei>\ ~ 1, the polynomial-type terms should 
resum, but it is hard (if not impossible) to predict the resummed form. We will therefore ignore these terms and 
concentrate on the exponential, with the post hoc justification that the approximations appear to work very well at 
intermediate masses and redshifts where non-Gaussian A^-body simulations have been performed. 
We are then after the quantity 



/ng,oxc(j^, £i, £2, • ■ .) 



— exp 

TT 



TT 



(-1)" 



n—3 



dX 



2tt 



(A3) 



where we defined 



(f){X) = iXv + y^ j — e„_2 ; eo = 1 . 



(A4) 



The integral in Eqn. (jA3[) can in principle be done using a saddle point approximation with calculable corrections, 
e.g. along the lines presented in Dll. Unfortunately, as Dll showed, a perturbative treatment will necessarily require 
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the condition \€v\ < 1, which is not sufficient for our purpose. Going beyond this technical barrier requires knowledge 
of an infinite number of connected moments of the density field. A key simplification occurs if we assume (inspired 
by the heirarchy ej ^ (/nlV^)"') the exact relations 



ei , 



(A5) 



One could treat this as defining a specific model of non-Gaussianity. Throughout, we will numerically compute £i as 
in the local model of NG. The series in 0(A) can then be re-arranged to bring the function into closed form, 



1 



(/)(A) ^iXi^ + — {e-'^" +iXe-l 



with derivatives 



7(> 



0"(A) 



-iXe 



A saddle point now exists at A = A, where ^'(A*) = 0, i.e. 



= 1 



(A6) 



(A7) 



(A8) 



provided we have (p"{X^,) = —(1 + eu) < 0, i.e. if ei/ > —1. The leading saddle point result for the multiplicity follows 
from setting / dX/V^e'^^^^ = e-^^^'^ (|0"(A*)|)"i/2^ which gives 



/i 



NG,cxc 



— (1 + eiy) ^^'^ exp 

/NG,app(i^, e) , 



4 (ei' - (1 + ei^) ln(l + ei^)) 



(A9) 



the subscript "app" reminding us of the approximations involved in the derivation. 

One can also estimate the error involved in the saddle point approximation. The calculation proceeds along the 
lines discussed in Appendix D of Dll, by first using the series representation of 0(A) in the integral of e'^^^\ and then 
Taylor expanding to get (after a change of variables) the exact result 



dX 

, ( 

V2^ 



(1 



-1/2 



dy 
, — ( 



3! 



rl/2. 



4! 



1 / {-iyf 



3! 



T + .. 



(AlO) 



where r = e^/(l + ev) > 0. Clearly, terms involving half-integer powers of r will not contribute, and hence the relative 
correction to Eqn. ()A9|) is of order 0{t) and is calculable in principle to arbitrary order in r. Since e ~ 3 x 10~^/nl, 
the approximation remains valid even for considerably large positive values of /nl, and becomes increasingly better 
for large v. Negative /nl values do not fare as well though, due to the restriction ei' > —1. See the main text however 
for a simple way of extending the mass function to /nl < 0, which appears to work well in the regime where A^-body 
simulations have been performed. 

It is straightforward to check that for fixed /nl, the multiplicity /NG,app(j^, e) is monotonically decreasing with v 
at large v, and the mass function is therefore stable in its extreme tail, as needed. One can also easily check (using 
ln(l + x) = — X^j^il^^)"/"-) that the limit e — > recovers the Gaussian Press-Schechter result, 



lim/NG,app(i^, e) = \ -ve ^^"^ = /ps(i^) . 
In fact, retaining the leading terms in e recovers the MVJ exponential prefactor 



/NG,app(!^, e) 



(All) 



(A12) 



Since the Press-Schechter result is known to perform badly when compared with A-body simulations, we will follow 
the usual practice in the literature and prescribe a ratio of non-Gaussian and Gaussian mass functions. Assuming 
that the Gaussian mass function is well described by the Tinker et al. form (29j , the final mass function to use is given 
by Eqn. ([3]), with a modified definition of v as discussed in the main text. 
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